Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN I^TeX style file v2.2) 



SCUBA polarisation observations of the magnetic fields in 
the prestellar cores L1498 and L1517B 

J. M. D. Ward-Thompson^'^ R. M. Cmtcher^ 

^School of Physics and Astronomy, Cardiff University, 5 The Parade, Cardiff, CF24 3YB 

^Department of Astronomy, University of Illinois, 1002 West Green Street, Urbana, Illinois, IL61801, USA 

^ Observatoire de Bordeaux, 2 rue de I'Observatoire, 33270 Floirac, France 

Accepted 2006 March ?; received 2006 March 24; in original form 2006 January 17. 

ABSTRACT 

We have mapped linearly polarized dust emission from the prestellar cores L1498 and 
L1517B with the James Clerk Maxwell Telescope (JCMT) using the Submillimetre 
Common User Bolometer Array (SCUBA) and its polarimeter (SCUBAPOL) at a 
wavelength of 850 /xm. We use these measurements to determine the plane-of-sky 
magnetic field orientation in the cores. In L1498 we see a magnetic field across the 
peak of the core that lies at an offset of ~19° ± 12° to the short axis of the core. 
This is similar to the offsets seen in previous observations of prestellar cores. To the 
southeast of the peak, in the filamentary tail of the core, we see that the magnetic field 
has rotated to lie almost parallel to the long axis of the filament. We hypothesise that 
the field in the core may have decoupled from the field in the filament that connects 
the core to the rest of the cloud. We use the Chandrasekhar-Fermi (CF) method to 
measure the plane-of-sky field strength in the core of L1498 to be ^10 ± 7 fiG. 

In L1517B we see a more gradual turn in the field direction from the northern 
part of the core to the south. This appears to follow a twist in the filament in which 
the core is buried, with the field staying at a roughly constant ~25° ± 6° offset to 
the short axis of the filament, consistent with previous observations of prestellar cores. 
Hence these two clouds in an apparently similar evolutionary state, that exhibit similar 
masses, morphologies and densities, have very different magnetic field configurations. 
We again use the CF method and calculate the magnetic field strength in L1517B to 
be '^30 ± 10 fiG. Both cores appear to be roughly virialised. Comparison with our 
previous work on somewhat denser cores shows that, for the denser cores, thermal 
and non-thermal (including magnetic) support are approximately equal, while for the 
lower density cores studied here, thermal support dominates. 
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1 INTRODUCTION 

The exact manner in which stars form is still a matter of 
debate. One important aspect of this debate centres around 
the morphology and magnitude of the magnetic field within 
molecular clouds and the magnitude of its effect on the star 
formation process. Observations of the magnetic field before 
protostellar collapse are necessary to address this question, 
ideally in a sample of molecular clouds that do not contain 
protostars but are on the verge of protostellar collapse. 

Myers and co-workers developed a catalogue of dense 
molecular cores within dark clo uds, which they observ ed in 
CO and later in a mmonia (see: iBenson fc Mverslll98gl and 
references therein). iBeichman et alTl^S^ cross-referenced 
the Myers ammonia catalogue with the IRAS point source 



catalogue and identified starless cor es as those that had not 
yet fo rmed an embedded protostar. IWard- Thompson et alJ 
il994h observed the submillimetre continuum emission from 
starless cores and identified a subset of starless cores called 
prestellar (formerly pre-protostellar) cores that are thought 
to be about to form a protostar. 

iKirk. Ward- Thompson fc Andrei J2005l hereafter 
KWA) observed a sample of 52 candidate prestellar cores 
with the Submillimetre Common User Bolometer Array 
(SCUBA) camera on the James Clerk Maxwell Telescope 
(JCMT). Twenty-nine of the cores were detected and were 
separated into 'bright' and 'intermediate' groups based on 
their 850-^m peak fiux densities. The bright cores were 
see n to be more centrally cond e nsed. Following the method 
of IJessop fc Ward-ThompsonI (1200 J) . KWA showed that 
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the estimated lifetimes of the two groups were consistent 
with them being part of an evolutionary sequence. 

Ward-Thompson et al. (2000) presented the first ob- 
servations of the magnetic field geometry from prestellar 
cores when they mapped L183, L1544 and L43 with the 
SCUBA polarimeter. Their maps showed relatively smooth 
and uniform magnetic fields over the central core regions 
that were at an angle of a pproximately 30 ° to the projected 
minor axis of the cores. ICrutcher et ahl l|20^) used the 
IWard- Thompson et alJ i2000h data to estimate the strength 
of the magnetic fi eld in the plan e of the sk y using an effect 
first suggested by IChan drasckh ar fc Fermi il953l hereafter 
CF). 

CF postulated that when the magnetic field is frozen 
to the matter, turbulent motions will cause local perturba- 
tions in a uniform magnetic field. These perturbations to the 
magnetic field should manifest themselves as a random scat- 
ter superimposed on an otherwise regularly oriented mag- 
netic field patern. If the magnitude of the random compo- 
nent of the dispersion, 5(j>, and the turbulent motion, SVnt, 
are known then CF showed that it is possible to infer the 
strength of the uniform plane-of-sky field Bpos in a medium 
of density p, such that: 
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where p = mn{H2), m is the mean particle mass, n{H2) is 
the volume number density and AVnt ~ SVnt V 8 In 2 is the 
full- width at half maximum (FWHM) of the non-thermal, or 
turbulent, linewidth. Q is a calibration facto r, fixed by com- 
pariso n with simulations, for which we follow lCrutcher et alJ 
||200J) in adopting a value of 0.5. 

A crucial parameter in the study of magnetically influ- 
enced star formation is M/$, which is the ratio of the mass 
to the magnetic flux and is normally expressed as A, the ratio 
of the observed mass-to-flux ratio (M/$)o(,3 to the critical 
mass-to-flux ratio {M/^)crit- Above the critical mass-to-flux 
ratio the magnetic field cannot support the observed mass 
against gravitational collapse. 

Inferring A from observations is possible if the column 
density and the magnetic field strength B can be deter- 
mined such that: 



A = 



B 



(2) 



where m = 2.8mH, allowing for 10% He by number, with 
N( H2) in cm~^ and B in p G. 

ICrxitcher et alJ ||2Qq3) estimated that the three cores 
they analysed had Bpos ~ 80 — 160/iC and A ~ 1.9 — 2.6. 
When corrected downwards by a statistical factor of three 
for geometrical bias this showed that the magnetic field 
strength was such as to render each core slightly magnet- 
ically subcritical. 

However, both Ward- Thompson et al. (2000) and 
ICrutcher et al.l i2004l) pointed out that the offset between 
the core short axis and the magnetic field direction was not 
consist ent with magnetically-regulat ed star formation mod- 
els (e.g. lCiolek fc MouschoviaJl998l and references therein). 
These models predict preferential collapse along the short 
axis of the core, yielding a field direction parallel to the 
short axis. This was explained as a projection effect in tri- 



axial cores (Basu 2000; Ciolek & Basu 2004), and it was 
claimed that the data were consistent with the models. 

All three of the prestellar cores previously mapped with 
the SCUBA polarimeter belong to the group of bright cores 
in the KWA survey. In this paper we present data for two 
more cores, L1498 and L1517B, which are somewhat less 
luminous. L1517B is a bright core and L1498 is in the inter- 
mediate group, although both lie near the borderline, and 
they are in fact quite similar. The goal is to examine whether 
the magnetic field morphology and/or strength change with 
submillimetre brightness, mass or evolutionary stage. 

We describe the observational parameters in the next 
Section and present our data in Section |21 In Section 2] we 
apply the CF technique to the new data and calculate vari- 
ous physical parameters of the cores. 



2 OBSERVATIONS 

Submillimetre continuu m observations at 85 pm were car- 
ried out using SCUBA llHolland et al.| ll999) on the JCMT 
on Mauna Kea, Hawaii, during the evenings of 1999 August 
13-15 between the hours 01:30 and 09:30 HST (11:30-19:30 
UT) and on the evenings of 2004 February 2-3 between the 

hours 17:30 and 01:30 HST (03:30-11:30 UT)^ 

SCUBA was used with the SCUBAPOL ('Greav es et all 

l2003h polarimeter, which uses a rotating half-waveplate and 
fixed analyser. The waveplate is stepped through sixteen 
positions (each offset from the last by 22.5°) and a Nyquist- 
sampled image (using a 16-point jiggle patte rn) is taken at 
each waveplate position iCreaves et alJl2003D . The observa- 
tions were carried out whilst chopping the secondary mirror 
120 arcsec in azimuth at 7 Hz and synchronously detecting 
the signal, thus rejecting sky emission. 

The integration time per point in the jiggle cycle was 
1 second, in each of the left and right telescope beams of 
the dual-beam chop. The total on-source integration time 
per complete cycle was 512 seconds. Two sources were ob- 
served, L1498 and L1517B, and the entire observing process 
described above was repeated 46 times for L1498 and 32 
times for L1517B. 

The instrumental polarisation (IP) of each bolometer 
was measured on the planets Mars and Uranus. This was 
subtracted from the data before calculating the true source 
polarisation. The mean IP was found to be 0.93 ± 0.27%. It 
was found necessary to offset the pointing between repeated 
observations so that holes left by the exclusion of bolome- 
ters with significantly higher than average noise statistics 
could be filled in. The submillimetre zenith opacity for at- 
mospheric extinction removal was determined by compar - 
ison with the 1.3-mm sky opacity ([Archibald et al.ll2002f) . 
The average 850-/im zenith opacity was 0.28 during the 1999 
observations and 0.19 during the 2004 observations, corre- 
sponding to a mean zenith transmission of 76 and 83% re- 
spectively. 



3 RESULTS 

Figures Q and |21 show maps of the prestellar cores L1498 
and L1517B. The greyscale and contours show the Stokes I 
(intensity) of the dust emission. The absolute calibration of 
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the intensity is not maintained through polarimetry obser- 
vations. Therefore we simply show the contours as a frac- 
tion of the peak intensity and refer to previous observations 
to obtain the absolute intensity calibration and parameters 
derived therefrom. For example, KWA found peak flux den- 
sities of 120 and 170 mjy/beam and integrated flux densi- 
ties of 2.3 and 2.6 Jy respectively for L1498 and L1517B at S 
850 nm. We use these data below to calculate core masses g 
and densities. g 
We note that there is a 180° ambiguity in the field di- ]j 
rection, as two anti-parallel fields would generate the same 
polarisation. The correct mathematical term for such a vec- 
tor with a 180° ambiguity is a 'half vector'. Therefore, this 
is the term that we use. The half vectors in Figures 1 & 
2 show the direction of the magnetic field in the plane 
of the sky. This is inferred from the assumption that the 
dust grains are magnetically aligned and that the magnetic 
field direction is orthogona l to the direction of polarisation 
jPavis fc GreensteinlllQSll) . Therefore the polarisation half 
vectors have been rotated by 90° to indicate the field direc- 
tion. 

The half vectors in Figures 1 & 2 are binned at 12- 
arcsec spacing to roughly match the JCMT beam FWHM 
of ~14 arcsec at 850 fim. Only those half vectors with a 
signal-to-noise ratio >2 a are plotted. This corresponds to a 
maximum position angle uncertainty of ±14°. For each core 
the highest detection was found to be just under 8a and the 
mean half vector signal-to-noise ratio was ~4-5 a. 

Note that, whenever we quote an error-bar on a po- 
larisation half vector position angle measurement, we are 
quoting the observational error rather than the dispersion 
observed from one half vector to another. This error can be 
calculated for a single half vector from its signal-to-noise ra- 
tio, based on its stability as a function of time throughout 
the observations. This can be cross-checked using the known 
sensitivity of the instrument, the observing conditions and 
the total integration time. The error-bar for a group of half 
vectors can then be calculated from each half vector's indi- 
vidual error-bar. 

3.1 L1498 

Figure 1 shows an intensity map of L1498 at 850 fim as a 
greyscale with contours overlaid. Superposed on this are a 
series of polarisation half vectors which have been rotated 
by 90° to illustrate the magnetic field direction. The length 
of each half vector is proportional to the percentage po- 
larisation observed. Looking first at the morphology of the 
core, we see an elongated core extending northwest to south- 
east, with a brighter peak at its northwestern end and a tail 
extending to the southeast. KWA also mapped L1498 at 
850 /im, but in photometric mapping mode only, and ob- 
served a slightly more rounded morphology. 

The intensity map presented in Figure Q is deeper and 
larger in extent than the previous map and shows that the 
chop direction chosen for the original map was almost par- 
allel to the long axis of the core (the orientation of the core 
on these scales was not known at the time of the original ob- 
servations). The chop function of the KWA data (see their 
section 5.4) caused an underestimate of the structure paral- 
lel to the major axis of the core and an offset in the apparent 
position angle of the core. It also led to an underestimate of 



_, , , 1 

. 1 20% 


1 1 1 1 1 1 1 r 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 

^^^^^^^^^ o -_ 


























, , 1 





58 56 54 52 4:10:50 48 46 



Right ascension (2000) 

Figure 1. The prestellar core L1498 seen in submillimetre contin- 
uum dust emission at a wavelength of 850 /^m. The contours are at 
40, 60, 80 & 90% of peak intensity. The polarisation half vectors 
have been rotated by 90° to show the magnetic field direction on 
the plane of the sky. The length of each half vector is proportional 
to the percentage polarisation. The ovals show groups of half vec- 
tors with roughly similar orientations and the line through each 
oval shows the approximate mean weighted PA of each group. The 
crosses mark th e centr es of integrated CCS emission observed by 
iLai fc CrutcheJ i2000D . The star marks the location o f the CCS 
Zeeman observation undertaken bv lLevin et al.l 1200 Jl . 



the magnitude of the extended flux density. Correcting for 
the chop function takes L1498 to being a borderline 'bright' 
core, even more similar to L1517B. 

The map of L1498 in Figure Q shows an elliptical core 
with a position angle (PA) for the long axis of the extended 
emission of ~135° ± 5°. The flux distribution is asymmetric 
along the long axis as mentioned above. The far-infrared 
maps made by the Infrared Space Observatory (ISO) of this 
region show that the L1498 core is at the northwestern end 
of a filament extending from the main cloud which lies to 
the southeast (Ward- Thompson, Andre & Kirk 2002) . 

The surface of this filament appears to be being 
heated by the inter-stellar radiation field (ISRF) and has 
a higher colour tem perature than the interior of the core 
JWard- Thompson et al. 2002) . The PA of the major axis of 
the core seen in our new data is coincident with the PA of the 
filament seen in the ISO maps (c.f. figure 7 of KWA). The 
tail of emission seen in Figure 1 heading to the southeast of 
the core is the start of this filament. 

However, there is an anti-correlation between the peak 
of the 850-/im emission and the peak of integrated CCS 
and CS emission. The latter peaks towards the SE of the 
core ( Tafalla ct al. 2004; Lai & Crutchcr 2000). Figure 8 of 
iLai fc Crutcheil 1I2OOOII shows what appears to be a CCS ring 
shape with a central depleted 'hole'. The peaks of the CCS 
emission are shown as crosses in Figure Q 

We detect twenty polarisation half vectors across L1498 
above the 2(t cut-off. There appears to be a dramatic change 
in field direction from the half vectors in the brighter north- 
western part of the core to those in the southeastern part. 
The two groups are shown encircled on Figure In between 
there is a gap in the half vectors where presumably the two 
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different field directions compete and cancel the net polari- 
sation. 

The thirteen half vectors on the core peak have a 
weighted mean PA of 64° ± 7°, which is at an angle of 19° ± 
12° to the PA of the short axis of the core. This alignment 
offset of the magnetic field to the orientation of the core 
short axis is similar to the ~20-30 ° offset found in previ- 
ous o bservations of prestellar cores (IWard-Thompson et alJ 
l2nn(i) . as discussed in section 1 above. 

The seven half vectors away from the core peak in the 
southeastern tail have a weighted mean PA of 124° ± 6° and 
are approximately parallel (11° ± 11° offset) to the long axis 
of the core. This rotation of the magnetic field through 60° 
± 11° is very marked and is most u nusual. It has not been 
obser ved in a prestellar core before (|Ward-Thompson et alJ 
l200[i^ . We note that the seven half vectors in the core's tail 
are coincident with the 'hole' in CCS emission discussed 
above, and with a region of anomalous velocity seen in N2H"'" 
(Tafalla et al. 2004). 

However, we noted above that the tail to the southeast 
seen in Figure 1 is also seen in the ISO data to extend even 
further to the southeast, as a filament linking up to the rest 
of the cloud (see figure 7 of KWA). So it may be that the 
extended filament has a different B-field direction from the 
dense core at its head and that the two fields have essentially 
become decoupled. If that is the case, then it seems quite 
remarkable that the two fields have become so decoupled as 
to lie at such a large apparent angle to one another. 



3.2 L1517B 

The map of L1517B in Figure |5| shows a peak of emission 
with a tail that extends to the southwest. The composite 
ISO and SCUBA maps of KWA show that L1517B extends 
in a truncated filament northeast of the main L1517 cloud. 
The peak of submillimetre intensity is near the 'head' of this 
filament and the tail points back towards the main cloud. 
Hence it has a somewhat similar morphology and situation 
to L1498, with a core near the head of a filament. 

However, there is also some evidence in the intensity 
map of Figure 2 that there is a slight extension of the fila- 
ment to the north of the peak, and that the filament turns 
slightly. Furthermore, the ISO data also show some evidence 
for a twist in the filament just to the north of L1517B (see 
figure 7 of KWA). Hence we see that the elongated filament 
in which the peak sits turns slightly from its northern end to 
its southern part. This is illustrated in Figure 2 by a curved 
dashed line which follows the ridge of the filament. 

We detected seventeen polarisation half vectors across 
L1517B above the 2-a cut-off. They have a weighted mean 
PA of 106° ± 7°. However, the morphology of the half vec- 
tors does not at first appear to be a simple uniform pattern. 
They were split into three groups of similar orientation half 
vectors - the northern four half vectors (PA ~82° ±6°), the 
middle nine half vectors (PA ~94° ±6°), and the southern 
four half vectors (PA ^156° ± 5°) - as shown in Figure |5| 
There appears to be a progression in mean PA from the 
north to the south, indicating a change in the mean field 
direction. 

The PA of the long axis of the filament was measured at 
the centre of each oval, and found to be 15° ± 5°, 35° ± 5° 
and 50° ± 5° from north to south respectively. Therefore the 
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Figure 2. The prestellar core L1517B in submillimetre dust emis- 
sion at 850 fim. The contours are at 40, 60, 80 & 90% of peak 
intensity. The polarisation half vectors have again been rotated by 
90° to show the plane-of-sky magnetic field direction. The ovals 
show groups of half vectors of similar orientation and the lines 
through the ovals show the approximate mean weighted PA of 
those half vectors. The dashed line shows the ridge of the larger 
filament in which the core is embedded. 



angle between the field direction and the filament elongation 
axis is 67° ± 11°, 59° ± 11° and 75° ± 10° from north to 
south respectively. 

The mean of these is 65° ± 6°, and hence the mean of 
the offset from the short axis is 25° ± 6°, consistent with 
that seen in the core of L1498, and with previous observa- 
tions of prestellar cores ( Ward-Tliomoson ct al. 2000). So we 
see that the field is simply following the twist in the filament 
and remaining at a roughly constant offset to it, as was seen 
before. However, we note that the sign of the offset flips in 
the south. This flip cannot be explained by the argument 
that this is simply a projection effect in a triaxial core, as 
suggested by the proponents of the magnetically regulated 
models (Basu 2000; Ciolek & Basu 2004). 



4 DISCUSSION 

We can use our measurements not only to map the orienta- 
tion of the magnetic field, but also to estimate the strength 
of the field in the plane of the sky, using the CF method, as 
discussed in section 1 above. For this we need to measure the 
dispersion in the polarisation half vectors (corrected for dis- 
persion caused by observational errors), and to compare this 
with the mass, density and non-thermal spectral linewidth 
in each core. 

The densities of the cores were estimated by assuming 
that they are triaxial ellipsoids, with the line-of-sight axis 
being the mean of the other two axes. We use the continuum 
data of KWA to calculate the dust mass, and hence the total 
mass of each core. The conversion from total mass to mass of 
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Source 


L1498 


L1517B 


Distance (pc) 


140 


140 


Diameter (pc) 


0.04 X 0.04 


0.03 X 0.02 


• OOS \ / 


0.3 ± 0.1 


0.3 ± 0.1 


N{H2) (10^2 cm-2) 


1 ± 0.5 


3 ± 2 


n{H2) (10^ cm~3 ) 


1 ± 0.5 


6 ± 3 


5^ n 


40 ± 7 


30 ± 5 


Bpos (fJ-G) 


10 ± 7 


30 ± 10 


^oba 


8 ± 5 


7 ± 4 


^corr 


3 ± 2 


2 ± 1 


MB,crit (Mq) 

Ay (kms-1) 


0.1 ± 0.07 


0.15 ± 0.1 


0.204 ± 0.01 


0.195 ± 0.01 


MaV (Mq) 


0.8 ± 0.13 


0.4 ± 0.15 


Afs,„ir (Mq) 


0.9 ± 0.2 


0.55 ± 0.25 


A/p,ob3 (Mq) 


0.8 ± 0.4 


0.6 ± 0.3 



Table 1. Calculated parameters of the prestellar cores L1498 and 
L1517B. For the distances we followed KWA. The sizes were mea- 
sured from Figures 1 & 2. The masses M^jj and column densities 
N{H2) were estimated from the submillimetre photometric map- 
ping of KWA. The volume densities n{H2) were calculated assum- 
ing that the line-of-sight dimension of each core was the mean of 
its other two dimensions. The half vector dispersion 5<f> was cor- 
rected for the dispersion due to observational error. The plane-of- 
sky B-field strength Epos was calculated from these parameters 
by the CF method. The ratio of the observed mass-to- flux ratio to 
the critical mass-to-flux ratio, Ao;,^, was calculated from the mass 
and magnetic field strength, and corrected by a mean geometrical 
factor of 3 to produce the corrected value \corr- The critical mass 
that could be supported by the field Ms^crit was calculated from 
this. The linewidths Ay are from lTafalla e'tli l. 1 2004). The (ther- 
mal plus nonthermal) virial mass M^v was calculated from the 
linewidth. The total virial mass Mg was calculated by adding 
together Mg ^rit and M^v- both cases Mg ^i^ is larger than 
Mobs} although M^bs must be corrected for external pressure to 
obtain Mp ^bst which is seen to be consistent with Mg^yir in both 
cases (see text for discussion). 

H2 is made using the standard solar abundance of H. These 
values were calculated and are listed in Table 1. 

The half vector dispersion S(j) was corrected for the dis- 
persion due to observational error (see: ICrutcher et al.l2004l 
for discussion). Bpoa was calculated from these parameters 
by the CF method. The ratio of the observed mass-to-flux 
ratio to the critical mass-to-flux ratio, A, was calculated 
from the mass and magnetic field strength , and corrected 
by a m ean geometrical factor of 3, following ICrutcher et alJ 
J2004I) . 

For L1498 we calculated the size, mass and density of 
the core peak only. To do this we measured the flux density 
in the northern circular aperture on Figure 1 using the KWA 
data. We corrected for the chop function along the long axis 
of the core, as discussed in section 3.1 above. We used this to 
calculate the column density and hence the mass. Similarly 
we only measured the half vector dispersion of the thirteen 
half vectors on the core peak, so as to exclude from our cal- 
culations the large turn of the apparent field direction from 
the core peak to the filamentary tail discussed in section 3.1 
above. 

Our calculations using the CF method give us a plane- 
of-sky magnetic field strength for L1498 of ~10 ± 7 /iG. 
This yields a geometrically-corrected ratio, Xcorr, of the 
observed to critical mass-to-fiux ratios of 3 ± 2. Zeeman 



measur ements of the line -of-sight magnetic field strength in 
L1498 iLevin et alJl200ll ) were taken at the central position 
of the core, marked by a star on Figure 1. These gave a field 
strength of 48 ± 31 nG in the line of sight. These two num- 
bers are comparable, to within errors, so we are not seeing 
the field at any particularly special angle. Adding the two 
components in quadrature gives a total field strength of 49 
± 33 /iC Hence we see that our geometrical correction fac- 
tor of ~3 is consistent to within the errors, and we therefore 
leave it unchanged to preserve comparability between the 
results for different cores. From our derived value of Xcorr 
we calculate the magnetic critical mass to be 0.1 ± 0.07 M0. 

For L1517B we calculated the half vector dispersion sep- 
arately in each of the three ovals shown in Figure 2 above, 
and then took the weighted mean of the three dispersions. 
Thus we excluded any contribution to the dispersion from 
the large-scale twist in the magnetic field discussed in sec- 
tion 3.2 above. Hence the CF method yields a plane-of-sky 
magnetic field strength for L1517B of ~30 ± 10 fiG. This 
gives a geometrically-corrected ratio, Xcorr, of the observed 
to critical mass-to-flux ratios of 2 ± 1. We note that in both 
cores the half vector dispersion is greater than the recom- 
mended upper limit of 25° for flnding accurate CF results 
(Ostriker, Gammie & Stone, 2001), so our values may be 
somewhat more uncertain as a result of this. We obtain a 
magneti c critical mass for L1517B of 0.15 ± 0.1 Mq. 

Taf afla et all 1I2OO4) observed L1498 and L1517B in 
multiple transitions of a range of nitrogen- and carbon-based 
molecules. They recorded respective mean intrinsic N2H''" 
FWHM linewidths of 0.204 and 0.195 kms~^ respectively, 
as listed in Table 1. We note that these authors found the 
same linewidths for NH3. They also found the non-thermal 
linewidths AVnt to be 0.12 and 0.11 kms~^ respectively, 
as used in the CF calculation above. These non-thermal 
linewidths are approximately a third to a half of thos e found 
for the three cores studied bv lCrutcher et aD ll2004l) . 

We used the total linewidths to calculate the (ther- 
mal plus non-thermal) virial masses of the cores, Mav , that 
could be supported virially by these linewidths, and found 
them to be 0.8 and 0.4 M© for L1498 and L1517B respec- 
tively. We list these in Table 1. We note that the mass able 
to be supported in this way is greater than the observed 
mass in each case. 

The parameters listed in Table 1 show that L1498 and 
L1517B are rather similar in terms of their properties such 
as size, mass, density and velocity dispersion. Therefore, the 
fact that their fleld morphologies are very different from one 
another, as discussed in section 3 above, is somewhat sur- 
prising. Table 1 shows that both cores appear to exceed their 
critical mass-to-fiux ratios by a factor of ~2-3. We also see 
that the (thermal plus non-thermal) virial mass, Mav, in 
each case is larger than the observed mass. 

We can estimate the total virial mass, M_g,t,ir, (exclud- 
ing the effects of external pressure) by adding together the 
(thermal plus non-thermal) virial mass, Mav, and the mag- 
netic critical meiss, Mg^crit- We find values of 0.9 and 0.55 
M0 respectively. We list these in Table 1. In both cases we 
see that the total virial mass, Ms.uir, is larger than the ob- 
served mass of the core. 

In the foregoing we have simply calculated the masses 
in the regions for which we have measured the magnetic field 
strengths. However, we know from KWA that the total mass 
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of each core is greater than simply the mass of the central 
region. Hence we account for the extra mass by considering 
this as an external pressure acting on the central region in 
each case. 

Wc can estimate this external pressure for a spherically 
symmetric core using the expression AwR^nextkT, where k 
is Boltzmanns's constant, T is the core temperature (found 
from KWA), riext is the external surface density and R is 
the radius under consideration, next was calculated by mea- 
suring the difference in mass between apertures of radii R 
and R + dr and assuming that this mass came from a shell 
of material of thickness dr. The density in this shell was 
found to be ~ 1 X lO'"' cm"-^ for L1498 and - 3 x lO''' cm"-^ 
for L1517B. The value for L1498 is in fact consistent with a 
uniform density out to the maximum radius probed. 

We can incorporate this, as a correction for external 
pressure, to our observed mass Mobs, to obtain Mp^ois (c.f. 
Ward- Thompson 2002; Ward- Thompson ot al. 2006). We 
find that the additional pressure translates into an equiv- 
alent additional mass of 0.5 and 0.3 Mq in L1498 and 
L1517B respectively. Adding these to Mobs we obtain values 
of Mp,obs of 0.8 and 0.6 M© in L1498 and L1517B respec- 
tively. These values arc listed in Table 1. We estimate that 
the error-bars on these values could be as high as ~ 50%, 
nevertheless we see that MB,vir ~ Mp^obs for both cores, and 
thus both are consistent with being virialised. 

These two cores appear somewhat different from the 
three cores studied by Crutcher et al. (2004), in that they 
have signicantly lower non-thermal linewidths and mag- 
netic field strengths. We commented above that the three 
cores previously studied (L1544, L183 & L43) were all in 
the 'bright' core category of KWA, whereas L1498 and 
L1517B studied here both lie near the border-line between 
the 'bright' and 'intermediate' cores of KWA. Therefore we 
postulate that perhaps we are seeing a trend for the denser 
'bright' cores to have higher non-thermal linewidths and 
magnetic field strengths. In addition, for the denser cores 
thermal and non-thermal (including magnetic) support are 
approximately equal, while for the lower density cores, ther- 
mal support dominates. 



5 CONCLUSIONS 

We have mapped linearly polarized dust emission from 
the prestellar cores L1498 and L1517B at a wavelength of 
850 /im and used the measurements to determine the plane- 
of-sky magnetic field orientation in the cores. In L1498 we 
saw a magnetic field on the core peak that is offset from the 
core short axis by ~19° ± 12°. To the southeast of the peak, 
in the filamentary tail, we saw that the magnetic field has 
rotated to lie almost parallel to the long axis of the filament. 
We hypothesised that the field in the core may have decou- 
pled from the field in the filament that connects the core to 
the rest of the cloud. We used the CF method to measure 
the magnetic field strength in the core of L1498 and found 
it to be ~10 ± 7 /nG. 

In L1517B we saw a more gradual turn in the field di- 
rection from the northern part of the core to the south. This 
appears to follow a twist in the filament in which the core is 
buried, with the field staying at a roughly constant offset of 
~25° ± 6° to the short axis of the filament. The CF method 



was used to calculate the magnetic field strength in L1517B, 
which was found to be ~30 ± 10 fiG. Both cores were seen 
to be roughly virialised. 
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